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We study the roles of the tensor and pairing correlations on the halo formation in 11 Li with an 
extended 9 Li+n+n model. We first solve the ground state of 9 Li in the shell model basis by taking 
2p-2h states using the Gaussian functions with variational size parameters to take into account the 
tensor correlation fully. In n Li, the tensor and pairing correlations in 9 Li are Pauli-blocked by 
additional two neutrons, which work coherently to make the configurations containing the Opi/2 
state pushed up and close to those containing the lSi/2 state. Hence, the pairing interaction works 
efficiently to mix the two configurations by equal amount and develop the halo structure in 11 Li. 
For 10 Li, the inversion phenomenon of s- and p-states is reproduced in the same framework. Our 
model furthermore explains the recently observed Coulomb breakup strength and charge radius for 
"Li. 

PACS numbers: 21.60.Gx, 21.10.Pc, 27.20. +n 



I. INTRODUCTION 

A pioneering secondary-beam experiment found that 
the size of 11 Li was surprisingly large, which was outside 
of the common sense of nuclear physics [l|. This large 
size was later interpreted as due to the halo structure of 
two neutrons around the 9 Li core nucleus @- This find- 
ing together with others motivated the nuclear physics 
community to start a new research field for the study of 
unstable nuclei and built new facilities of Radioactive Ion 
Beams (RIB) in several laboratories as RIKEN, MSU, 
GSI, GANIL and others. Many experimental findings 
were shown later for 11 Li: a) The halo neutrons have al- 
most equal amount of the s-wave component with respect 
to the p-wave component Q . b) The dipole strength dis- 
tribution has a large enhancement near the threshold [|| . 
c) The charge radius is larger than that of 9 Li [1, @. 

The biggest puzzle from the theory side is the large 
s-wave component for the halo neutrons. If we interpret 
this fact in the shell model, the shell gap at N — 8 has 
to disappear. However, the mean field treatment of a 
central force is not able to provide the disappearance of 
the N = 8 shell gap. So far, there were many theoretical 

studies for n Li@, i, 1, M [LI E3, M M M EE El 

and essentially all the theoretical works of 11 Li had to ac- 
cept that the lsi/ 2 single particle state is brought down 
to the degenerated energy position with the Opi/2 state 
without knowing its reason[9(. It is therefore the real 
challenge for theoretician to understand this disappear- 
ance of the N = 8 shell gap, to be called the s-p shell 
gap problem, which is worked out in this paper by de- 
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veloping a framework to treat the tensor force explicitly 
in the nucleon-nucleon interaction. The halo structure of 
n Li is also related with the ls-state and the Op-state in 
10 Li. Several experiments suggest that the dual states of 
the Si/2-state coupled to the 3/2~ proton state appear 
close to the threshold energy of 9 Li+n together with the 
dual states of the pi/ 2 -state [lj| H(| [H|- This property 
is known as the inversion problem of s- and p-states seen 
in the N — 7 isotone[20j. 

The pairing correlation was considered important for 
the s-p shell gap problem 0, [l3j]. The additional neu- 
trons would act to block the pairing correlation of the 
core when one neutron in 10 Li or two neutrons in n Li 
are added in the p-orbit[H, HH, [23|. Our calculations 
with this pairing-blocking effect improve somewhat the 
descriptions of 10 Li and 11 Li, but it was not sufficient to 
explain the large s-wave component in n Li[l3|. We also 
pointed out the different roles of the pairing-blocking be- 
tween 10 Li and n Li. We need further mechanism for the 
increase of the s-wave component in n Li. 

The tensor force, on the other hand, plays an impor- 
tant role in the nuclear structure. For example, the con- 
tribution of the tensor force in the binding of 4 He is com- 
parable to that of the central force 0, [25[. The tensor 
correlation induced by the tensor force was demonstrated 
to be important for the 4 He+n system (2|| H3, HH ■ In our 
recent study [29| . we developed a theoretical framework 
of the tensor-optimized shell model to treat the tensor 
force in the shell model basis explicitly including Ip-lh 
excitations. We found that the (0s 1 / 2 )~ 2 (0pi/ 2 ) 2 exci- 
tation of proton-neutron pair has a special importance 
in describing the tensor correlation in 4 He (2|| H^, [3(| • 
In the 4 He+n system, because this 2p-2h excitation re- 
ceives the strong Pauli-blocking from the last neutron 
occupying the pi/2-orbit, a considerable amount of the 
P1/2-P3/2 splitting energy in 5 He is reproduced This 
Pauli-blocking effect from the p!/ 2 -orbit caused by the 
tensor force should be present also for 11 Li. 
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Hence, it is very interesting to study the effect of the 
tensor correlation together with the pairing correlation 
for the s-p shell gap problem in 11 Li. This is the purpose 
of this paper. To this end, we shall perform the configu- 
ration mixing based on the shell model framework for 9 Li 
to describe the tensor and pairing correlations explicitly. 
In particular, we pay attention to the special features of 
the tensor correlation. For 11 Li, we shall solve the con- 
figuration mixing of the 9 Li+n+n problem which treats 
both correlations, and investigate further the Coulomb 
breakup strength of 11 Li and other observables to see 
the effect of these correlations. We also investigate the 
inversion phenomena of 10 Li considering the tensor and 
pairing correlations. 



II. MODEL 



Coupled 9 Li+n+n model of n Li 



We shall begin with the introduction of the model for 
'Li, whose Hamiltonian is given as 



i?( 9 Li) = j>-t c +£\ 



(i) 



z<3 



Here, ti, to, and Vij are the kinetic energy of each nu- 
cleon, the center-of-mass (cm.) term and the two-body 
NN interaction consisting of central, spin-orbit, tensor 
and Coulomb terms, respectively. The wave function 
of 9 Li(3/2~) is described in the tensor-optimized shell 
model[28|, [29|. We express 9 Li by a multi-configuration, 



TV 



5 Li) = a i $ 



3/2" 



(2) 



where we consider up to the 2p-2h excitations within the 

Op shell for <P/ in a shell model type wave function, 
and N is the configuration number. Based on the pre- 
vious study of the tensor-optimized shell model [28l. |29|. 
we adopt the spatially modified harmonic oscillator wave 
function (Gaussian function) as a single particle orbit 
and treat the length parameters b a of every orbit a of 
Os, Opi/2 and Ops/2 as variational parameters. This vari- 
ation is shown to be important to optimize the tensor 
correlation [H [H M, EH . 

Following the procedure of the tensor-optimized shell 
model, we solve the variational equation for the Hamilto- 
nian of 9 Li and determine {c^} in Eq. ^ and the length 
parameters {b a } of three orbits. The variation of the 
energy expectation value with respect to the total wave 
function \&( 9 Li) is given by 



(#|J3"( 9 Li)|*> 







which leads to the following equations: 
5(*|H( 9 Li) - E\V) <9(*|#( 9 Li 



db a 



0, 



E\9) 



da; 



(3) 



0. (4) 



Here, E is the total energy of 9 Li. The parameters {b a } 
for the Gaussian bases appear in non-linear forms in the 
total energy E. We solve two kinds of variational equa- 
tions in the following steps. First, fixing all the length 
parameters b a , we solve the linear equation for {a.;} as 
an eigenvalue problem for H( 9 Li). We thereby obtain 
the eigenvalue E, which is a function of {b a }. Next, we 
try various sets of the length parameters {b a } to find 
the solution which minimizes the energy of 9 Li. In this 
wave function, we can optimize the radial form of single 
particle orbit appropriately so as to describe the spatial 
shrinkage of the particle state, which is important for the 
tensor correlation [H, [U, HO, S3 • 

For 11 Li and 10 Li, their Hamiltonians are written in 
terms of 9 Li+n+n and 9 Li+n, respectively, and are given 
as 



2 2 

(3) 



k=0 



k=l 



(5) 



tf( 10 Li) = i/( 9 Li) + ]TT fe -4 2) + V/ c „, (6) 



k=Q 



where H( 9 U), T k , and T^ 2) are the internal Hamil- 
tonian of 9 Li given by Eq. ([1]), the kinetic energies of each 
cluster (k — for 9 Li) and the cm. terms of three or two 
cluster systems, respectively. V cn ^ are the 9 Li core-n in- 
teraction (k — 1,2) and V nn is the interaction between 
last two neutrons. The wave functions of 11 Li and 10 Li 
with the spin J and J', respectively, are given as 



N 

vF( n Li) = Yl A {^i /2 ~>^°( nn ^}> ( ? ) 

i 

N 

'( 10 Li) - £^{[^xf>)] J '}. (8) 



We obtain the coupled differential equations for the neu- 
tron wave functions x Jo ( nn ) an( i X J °{ n )i where Jo and 
Jg are the spins of the additional neutron part of 11 Li 
and 10 Li, respectively. To obtain the total wave function 
fy J ( Li) and *£ J ( 10 Li), we actually use the orthogonal- 
ity condition model (OCM) [7|, [l3|, [32j to treat the an- 
tisymmetrization between last neutrons and neutrons in 
9 Li. In OCM, the neutron wave functions x are imposed 
to be orthogonal to the occupied orbits by neutrons in 

9 Li, which depend on the configuration <E>^ 2 in Eq. ((2]). 
We obtain the following coupled Schrodinger equations 
with OCM for the set of the wave functions {Xi°( nn )} 



for Li and {Xi°( n )} f° r Li, where i — 1 



N 

E 



( 9 Li) 



k=l 



(9) 
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N 

E [( T ™i + v ™ + A >) s a + M 9Li ) 

3=1 

x x f°(n) = E X f>(n), (10) 
Ai = A £ MfeXM (11) 

aS*i( 9 Li) 

where ^( 9 Li) = ($f/ 2 ~|H( 9 Li)|$J /2 ">. and T r (2) 

are the total kinetic energies consisting of the relative 
motions for n Li and 10 Li, respectively. A< is the pro- 
jection operator to remove the Pauli forbidden states tp a 
from the relative wave functions [2^. |33| . where i/j a is the 
occupied single particle wave function of the orbit a in 
9 Li. This Ai depends on the neutron occupied orbits in 

the configuration <E>^ 2 of 9 Li and plays an essential role 
to produce the Pauli-blocking in n Li and 10 Li, explained 
later. The value of A is taken large as 10 6 MeV in the 
present calculation in order to project out the compo- 
nents of the Pauli forbidden states into an unphysical 
energy region. Here, we keep the length parameters {b a } 
of the single particle wave functions as those obtained for 
9 Li. 

We explain the method of treating the orthogonality 
condition inclu ding the particle-hole excitations of 9 Li in 
more detail[H, [28|. When the neutron orbit in 9 Li is 
fully occupied, the orthogonality condition for the last 
neutrons to this orbit is given by Aj in Eqs. ^ and (fTU)) . 
When neutron orbits in 9 Li are partially occupied, such 
as in the 2p-2h states, the last neutrons can occupy these 
orbits with particular probabilities, which are determined 
by the fractional parentage coefficients of the total wave 
functions of 10,11 Li consisting of 9 Li and the last neu- 
trons. 

We describe the two neutron wave functions x i n 
Eq. for 11 Li precisely in a few-body approach of the 
hybrid-TV model 0, El ; 

Xi°{nn) = x{°{nn,£ v ) + xi {nn,i T ), (12) 

where £ v and £ T are V-type and T-type coordinate sets 
of the three-body system, respectively. The radial part 
of the relative wave function is expanded with a finite 
number of Gaussian basis functions centered at the ori- 
gin. We use at most 15 Gaussian basis functions with the 
maximum range parameter 30 fm to describe the loosely 
bound wave function of neutron halo[32j. 

Here, we discuss the coupling between 9 Li and the last 
neutrons, whose details were already explained in the 
pairing-blocking case[H, HH, |35| . We consider the case 
of 11 Li. In the present three-body model, the Pauli for- 
bidden states for the relative motion provides the Pauli- 
blocking effect caused by the last two neutrons [H, [23| . 
This blocking depends on the relative distance between 
Li and the two neutrons, and change the structure of 
9 Li inside Li, which is determined variationally to min- 
imize the energy of the n Li ground state. Asymptot- 
ically, when the last two neutrons are far away from 



9 Li (£ VT — > oo), the effects of antisymmetrization and 
the interaction between 9 Li and two neutrons vanish in 
Eq. Therefore, any coupling between 9 Li and two 
neutrons disappears and 9 Li becomes its ground state. 
Namely, the mixing coefficients {ai} are the same as those 
obtained in Eq. © : 

<f J ( n Li) ► [*( 9 Li), X ,a inn)] " , (13) 

5 V)T — «» N 

* »'!>;'-' • (14) 

i 

Therefore, it is easy to obtain the following asymptotic 
forms of {Xi° ( nn )} from the above two relations: 

Xi°(nn) ► a t - X Jo (nn), (15) 

where i = 1, ■ ■ ■ , N. Eq. (fT5"|) implies that the asymptotic 
wave function of two neutrons Xi° (nn) is decomposed 
into the internal amplitude a; of 9 Li and the relative 
wave function x Jo (nn). Eqs. (fT3|) - (fT5|) give the bound- 
ary condition of the present coupled three-body model of 
11 Li. Contrastingly, when the two neutrons are close to 
9 Li, the two neutrons dynamically couple to the config- 
uration < i ) ^ 2 of 9 Li satisfying the Pauli principle. This 
coupling changes {ai} of 9 Li from those of the 9 Li ground 
state, and makes the tensor and pairing correlations to 
be different from those in the isolated case. For 10 Li, the 
similar coupling scheme is considered. The dynamical 
effect of the coupling arising from the Pauli-blocking is 
explained in the results in detail. 

B. Effective interactions 

We explain here the interactions employed in Hamil- 
tonians in Eqs. ((T|), ([5]) and ©. Before explaining the 
present interactions, we give a brief review of the sit- 
uation of the treatment of the effective interactions for 
the study of 9,10,11 Li. As was mentioned, most theo- 
retical studies based on the three-body model of 11 Li 
employ the state-dependent 9 Li-n potential where only 
the s-wave potential is made deeper than other partial 
waves 9], while the 9 Li core is described as inert. This 
state-dependence in the 9 Li-n potential is phenomeno- 
logically determined in order to satisfy the experimental 
observations of a large s 2 component and a two-neutron- 
separation energy of 11 Li, and a virtual s-state in 10 Li, 
simultaneously. On the other hand, for the nn part, the 
interaction having a mild short-range repulsion[lllll8j or 
the density-dependent one are often usedQ. However, 
even in the microscopic cluster models using an unique 
effective ^V^V interaction consisting of the central and LS 
forces [13 GJ], the s-p shell gap problem in n Li and 10 Li 
cannot be solved simultaneously. From these results, we 
consider that the usual approach based on the effective 
central and LS interactions may be insufficient to ex- 
plain the exotic structures of 10 ' n Li. For this problem, 
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even the so-called ab-initio calculations using the real- 
istic NN interactions, such as Green's function Monte 
Calro[36j, do not provide good results for n Li. 

In this study, we focus on the tensor correlation, which 
is newly considered to figure out the s-p shell gap prob- 
lem. To do this, we extend the three-body model of n Li 
to incorporate the tensor correlation fully, in particular, 
for the 9 Li part. In the present study, our policy for the 
study of n Li is to use the experimental informations and 
the corresponding theoretical knowledge for 9 Li and 10 Li 
as much as possible. Following this policy, we explain 
our interactions in three terms; Vij of H( 9 L\) in Eq. H}, 
core-n V cn and n-n V nn of the Hamiltonians in Eqs. (|5|) 
and ©. 

For the potential V nn between the last two neutrons, 
we take a realistic interaction AV8' in Eq. ©. Our in- 
terest is to see the n-n correlation in the two-neutron 
halo structure, and therefore it is necessary to solve two- 
neutron relative motion without any assumption. For 
this purpose, our model space of two neutrons using the 
hybrid -TV model shown in Eq. (|12[) has no restriction 
and wide enough to describe the short range correlation 
under the realistic nuclear interaction AV8'. Therefore, 
there is no parameter in the potential V nn . 

The 9 Li-n potential, V cn , in Eqs. © and ^ is 
given by folding an effective interaction, the MHN 
interaction]^, l3a |. which is obtained by the G-matrix 
calculation and frequently used in the cluster study of 
light nuclei @, H [H \M, II • In the 9 U+n system, 
the folding potential for the 9 Li density calculated by us- 
ing H.O. wave function has been discussed to reproduce 
the proper energies of the 10 Li spectra0, HI, H^. Fur- 
thermore, considering the small one-neutron-separation 
energy of 9 Li and a long-range exponential tail of the 
density, we improve the tail behavior of the folding po- 
tential to have a Yukawa type form [3, [3f|. Any state- 
dependence is not used in the present 9 Li-n potential, 
such as a deeper potential for the s-wave. This is possi- 
ble because the Pauli blocking effect of the single parti- 
cle state is in action and the state with the P1/2 orbit is 
pushed up in energy and becomes close to the state with 
the S1/2 statefTH, [23[. We will discuss the results on 10 Li 
after the discussion on 11 Li. We introduce one parameter, 
S, which is the second-range strength of the MHN poten- 
tial in the calculation of the 9 Li-n potential to describe 
the starting energy dependence dominantly coming from 
the tensor force in the G-matrix calculation [32, l38|. In 
the present calculation, we chose this d parameter to 
reproduce the two-neutron-separation energy of n Li as 
0.31 MeV after working out the tensor and the pairing 
correlation effects as explained later. It is found that this 
folding potential also reproduces the positions of the p- 
wave resonances in 10 Li, just above the 9 Li+n threshold 
energy [l9j. as shown in the results. 

Now we discuss the choice of the interaction between 
nucleons in the 9 Li core; Vij in 77 ( 9 Li), where we use the 
shell model wave functions for the 9 Li core in Eq. ^j. 
Since our main interest in this work is to investigate 



the role of the tensor force on the two-neutron halo for- 
mation, we describe the tensor correlation in addition 
to the pairing correlation in the 9 Li core based on the 
policy mentioned above. Along this line, recently we 
have many interesting works [30, Hi], [4(1 El, Wc 
have also studied the role of the tensor force in the shell 
model framework, and proposed the tensor-optimized 
shell model. [28|, l29( As a reliable effective interaction con- 
sidered from those studies, in this calculation, we use GA 
proposed by Akaishi [H, Ej, S3] for u y in Eqs. Q, © and 
©. This effective interaction GA has a term of the ten- 
sor force obtained from the G-matrix calculation using 
the AV8' realistic potential keeping the large momentum 
space Hi], In GA, the obtained 9 Li wave function in 
Eq. © shows smaller matter radius than the observed 
one due to the high momentum component produced by 
the tensor correlation[28|, [3(| [M|. Hence, we have to 
adjust the central force, which is done by changing the 
second range of the central force by reducing the strength 
by 21.5% and increasing the range by 0.185 fm to repro- 
duce the observed binding energy and the matter radius 
of 9 Li in the same manner as done for 4 He[28l[2^|. 

C. Tensor correlation in the Gaussian expansion 
method 

In the description of the tensor correlation, in princi- 
ple, we can work out a large space to include the full effect 
of the tensor force by taking 2p-2h states with very high 
angular momenta[2g|. In order to avoid large computa- 
tional efforts without loss of the physical importance in 
the result, we restrict the 2p-2h shell model states within 
the p-wave states for the description of 9 Li with a single 
Gaussian basis. We have studied that the superposition 
of the Gaussian bases improves the description of the spa- 
tial shrinkage for the particle states caused by the tensor 
correlation[28|,|29|. In this case, the wave function of the 
particle state ip a in 9 Li is expanded with a finite number 
of Gaussian basis functions in a jj coupling scheme as 

N a 

i> a = ^C Q ,„ CMa.n), (16) 
n=l 

Here n is an index for the Gaussian basis with the length 
parameter & Q .„. A basis number and the normalization 
factor for the basis are given by N a and Af a ,n, respec- 
tively. The coefficients {G QjTl } are determined variation- 
ally for the total wave function of 9 Li in Eq. ^j. Using 
this method, so-called the Gaussian expansion method 
(GEM) 43], the wave functions of the particle states are 
improved with an appropriate radial form, where the set 
of {b a , n } is suitably chosen [29ll43|. In particular, it was 
shown that the particle-hole excitations induced by the 
tensor force increase [29j. We have confirmed the GEM ef- 
fect on the (0s 1 / 2 ) _2 (0p 1 / 2 ) 2 component for 4 He in Fig.[TJ 
since the similar GEM effect is expected for 9 Li. As the 
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FIG. 1: Probabilities of the (Os 1/ / 2 )io (Opi/2)io component in 
4 He in the Gaussian expansion method are shown as a func- 
tion of the number of Gaussian basis. Two subscripts 10 rep- 
resents spin and isospin for the two-nucleon pair, respectively. 
The converged value is reproduced by enhancing the tensor 
matrix elements with one Gaussian basis by 50%. 
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FIG. 2: (Color online) Energy surface of Li with respect 
to the length parameters b a of Op orbits. The two minima 
indicated by (a) and (b) in the contour map correspond to the 
states due to the tensor correlation and the paring correlation, 
respectively. 



number of Gaussian basis increases for the particle states, 
the (0s!/2)~ 2 (0pi/2) 2 component increases and converges 
with three Gaussians. This converged value could be re- 
produced by increasing the matrix elements of the tensor 
force with a single Gaussian basis by 50% as shown in 

Fig.ru 

Similarly, the GEM effect also affects the Pauli- 
blocking caused by adding a neutron into the occupied 
neutron orbit in the core. In the scattering problem of the 
4 He+n system, we checked that the GEM effect on the 
Pauli-blocking is reproduced using the enhanced tensor 
matrix elements with a single Gaussian basis[28, 44, 45]. 
The Pauli-blocking could be considered to be almost pro- 
portional to the overlap between the wave functions of 
neutrons inside and outside the core. In this sense, the 
single particle properties of the particle states of the core 
can be described using the enhanced tensor matrix ele- 
ments. Therefore, in the present study, we adopt this 
enhanced tensor matrix elements with a single Gaussian 
basis in order to simulate the GEM effect. 



III. RESULTS 
A. 9 Li 

We first show the results of the 9 Li properties, which 
give a dynamical influence on the motion of last neutrons 
above the 9 Li core in n ' 10 Li. In Fig. ® we display the en- 
ergy surface of 9 Li as functions of the length parameters 
of two Op orbits, where &o s is already optimized as 1.45 
fm. There are two energy minima, (a) and (b), which 
have almost a common bo P3/2 value of 1.7-1.8 fm, and a 
small (0.85 fm) and a large (1.8 fm) bo Pl/2 values, respec- 
tively. The properties of two minima are listed in Table U 
with the dominant 2p-2h configurations and their prob- 
abilities. It is found that the minimum (a) shows a large 
tensor contribution, while the minimum (b) does not. 
Among the 2p-2h configurations, the largest probabili- 



TABLE I: Properties of 9 Li with configuration mixing. 
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Expt. 




(a) 


(b) 


(c) 




E [MeV] 


-43.8 


-37.3 


-45.3 


-45.3 


(V T ) [MeV] 


-22.6 


-1.8 


-20.7 




Rm [fm] 


2.30 


2.32 


2.31 


2.32±0.02[46] 


0p-0/i 


91.2 


60.1 


82.9 




(Op 3 /2)oi 2 (OPl/2)oi 


0.03 


37.1 


9.0 




(0Sl/2)ro 2 (0p 1 / 2 )?O 


8.2 


1.8 


7.2 





ties are given by (0s) 10 2 (0pi/ 2 ) 2 for (a), similar to the 

results in Ref. [H, [2!|, and (0p3/ 2 )oi 2 (0pi/2)oi; namely 
the Op shell pairing correlation for (b). These results in- 
dicate that the minima (a) and (b) represent the different 
correlations of the tensor and pairing characters, respec- 
tively. The spatial properties are also different from each 
other; the tensor correlation is optimized with spatially 
shrunk excited nucleons for (a) and the pairing correla- 
tion is optimized when two Op orbits make a large spa- 
tial overlap for (b) . In Table HI we show the results of 
the superposition of minima (a) and (b), named as (c), 
to obtain a 9 Li wave function including the tensor and 
pairing correlations, simultaneously. For (c), the favored 
two configurations in each minimum (a) and (b) are still 
mixed with the Op-Oh one, and the property of the tensor 
correlation is kept in (c). The superposed 9 Li wave func- 
tion possesses both the tensor and pairing correlations. 



B. Pauli-blocking effect in n Li 

We discuss here the Pauli-blocking effect in n Li and 
10 Li. We mainly consider the case of 11 Li, as shown in 
Fig. [3] For the 9 Li ground state (GS, upper panel), in 
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FIG. 3: Schematic illustration for the Pauli-blocking in n Li. 
Details are described in the text. 



TABLE II: 5 and the energy differences AE in MeV. 
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FIG. 4: (Color online) (Is) 2 probability P(s 2 ) and mat- 
ter radius R m of 1 Li with four models in comparison with 
the experiments ((a)H|, (b)[|(| and (c)H3]). The scale of 
P(s 2 ) (Rm) is right (left) hand side. 



addition to the Op-Oh state, 2p-2h states caused by the 
tensor and pairing correlations are strongly mixed. Let 
us add two neutrons more to 9 Li. When two neutrons oc- 
cupy the 0pi/2-orbit (middle panel), the 2p-2h excitations 
of the tensor and pairing correlations in 9 Li are Pauli- 
blocked, simultaneously [l3j. Accordingly, the correlation 
energy of 9 Li is partially lost inside n Li. For the (Is) 2 
case of two neutrons (lower panel), the Pauli-blocking 
does not occur and 9 Li gains its correlation energy fully 
by the configuration mixing with the 2p-2h excitations. 
Hence, the relative energy distance between (Op) 2 and 
(Is) 2 configurations of n Li is expected to become small 
to break the magicity in 11 Li. The same effect is also 
expected to explain the inversion phenomena of Is- and 
p-states in 10 Li. 

In order to confirm the above expectation of the block- 
ing effect on the (Is) 2 configuration of 11 Li, we dis- 
cuss here the configuration mixing including sd-shell 
for 9 Li. The (0p3/ 2 )oi 2 (ls)o 1 neutron pairing excita- 
tion in 9 Li is negligible 23], and the probability of the 
(0s)f 2 [(ls)(0d3/2)]io excitation induced by the tensor 
force is around 2%[29|]. The latter excitation is a proton- 
neutron pair, in which the ls-state is spatially shrunk 
about a half size of bo s due to the tensor correlation [29|. 
When the ls-state is occupied by a neutron in 9 Li, this 
ls-state brings a small overlap with the spatially ex- 
tended ls-orbit of the last neutrons in n Li. Then we 
have estimated that the blocking effect on the (Is) 2 com- 
ponent of 11 Li for this proton-neutron excitation is very 
small. We consider that the characteristics of the block- 
ing effect for 11 Li would not change, even if we include 
the sd-shell for 9 Li and n Li. 



C. n Li 

We perform the coupled three-body calculation of 
11 Li considering the tensor and pairing correlations fully, 
named as "Present" . In order to see the individual effects 
of the tensor and paring correlations, we also compare the 
results with other three kinds of calculations for 11 Li with 
different descriptions of 9 Li. "Inert core" is only the Op- 
Oh configuration of 9 Li. "Tensor" and "Pairing" are the 
ones in which the minimum (a) and (b) in Table U are 
adopted for 9 Li, respectively. For each calculation, we 
determine the parameter 8 in the 9 Li-n potential, shown 
in Table HD 

In Fig. QJ "Present" is found to give a large amount of 
the (Is) 2 probability P(s 2 ), 46.9% for the last two neu- 
trons and a large matter radius R m , 3.41 fm for Li, 
which are enough to explain the observations. The prob- 
abilities of (pi/ 2 ) 2 , (P3/2) 2 , {dz/2) 2 and (d 3/2 ) 2 for the 
last two neutrons are obtained as 42.7%, 2.5%, 4.1% and 
1.9%, respectively. In Fig. [4j when we individually con- 
sider the tensor and pairing correlations for 9 Li, P(s 2 ) 
is larger for the tensor case than for the pairing case. 
This means that the blocking effect from the tensor cor- 
relation is stronger than the pairing case. Finally, both 
blocking effects furthermore enhance P(s 2 ) and provide 
almost equal amount of (Is) 2 and (Op) 2 configurations. 
Hence, two correlations play important roles to break the 
magicity and make the halo structure for 11 Li. 

In Table HT1 we also estimate the relative energy differ- 
ence AE between (Is) 2 and (Op) 2 configurations for n Li 
using the mixing probabilities of these configurations and 
the coupling matrix element between them as 0.5 MeV 
obtained in Ref. [l3| . The present model is found to give 
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FIG. 5: (Color online) Calculated Coulomb breakup cross 
section measured from the 9 Li+n+n threshold energy. 

the degenerated energies enough to cause a large cou- 
pling between the (Op) 2 and (Is) 2 configurations by the 
pairing interaction between the last neutrons. 

In addition to the matter radius, the halo structure 
also affects the proton radius of 11 Li, because of the re- 
coil effect of the cm. motion. In the three-body model of 
11 Li, its proton radius (R p ) consisting of the proton ra- 
dius of 9 Li and the relative distance between 9 Li and the 
cm. of two neutrons (R c ^2n) with the following relation 

(i? 2 ( n Li)) = (i? 2 ( 9 Li)) + (Rl_ 2a ), (18) 

where the second term represents the recoil effect. When 
the halo structure develops, (i? 2 _ 2n ) 1S expected to be 
large. Experimentally, considering the nucleon radius, 
the charge radius of 11 Li was measured recently and its 
value is 2.467±0.037 fm, which is enhanced from the 
one of 9 Li, 2.217±0.035 fm @. The improved calcu- 
lation for the isotope shift determination Q shows that 
2.423±0.037 fm and 2.185±0.033 fm for n Li and 9 Li, 
respectively. The present wave functions provide 2.44 
fm and 2.23 fm for n Li and 9 Li, respectively, which 
are in good agreement with the experimental values. 
This enhancement is mainly caused by the large value 

°f \J (^c-2n) obtained as 5.69 fm. For comparison, the 
distance between last two neutrons is 7.33 fm, which is 
larger than the core-2n case. 

We further calculate the three-body Coulomb breakup 
strength of 11 Li into the 9 Li+n+ri system to investigate 
the properties of the dipole excited states and compare 
the strength with the new data from the RIKEN group [4J. 
We use the Green's function method combined with the 
complex scaling method [32] to calculate the three-body 
breakup strength[35j using the dipole strength and the 
equivalent photon method, where the experimental en- 
ergy resolution is taken into account [ij. We do not find 
any resonances with a sharp decay width enough to make 
a resonance structure. In Fig. [5J it is found that the 
present model well reproduces the experiment, in partic- 
ular, for low energy enhancement and its magnitude. See- 
ing more closely, however, our results seem to underesti- 
mate the cross section at E > 1 MeV, while overestimate 



at low energy peak region slightly. As a result, the inte- 
grated dipole strength for E < 3 MeV gives 1.35 e 2 fm 2 , 
which agrees with the experimental value of 1.42 ± 0.18 
e 2 fm 2 [|. 

For the reference, we calculate the strength with a po- 
tential model denoted as DR, in which the 9 Li core is 
inert and the 9 Li-n s-wave potential is deepened to re- 
produce 50% of P(s 2 ) in the n Li ground state. In this 
case, we obtain three dipole resonances of l/2 + , 3/2+ 
and 5/2+ states with 3/2~ ® 1", less than 0.5 MeV 
above the three-body threshold energy, similar to the re- 
sults of Ref. [Hj]. In our results, the 3/2+ state is lo- 
cated slightly lower than other two states, because of the 
J = 0~ = (L = 1) ® (S = 1) component for two neu- 
trons in Eq. ([7]), where L and S are the coupled angular 
momenta and spins of the last two neutrons, respectively. 
This component does not appear in the l/2 + and 5/2 + 
states. This difference for the dipole states makes a vis- 
ible splitting in the cross section before folding with ex- 
perimental resolution as shown in Fig. [5] If we fold the 
spectrum by the experimental resolution, two peaks are 
washed out but the strength distribution comes out dif- 
ferent from both the present result and the experiment. 
The detailed analysis of the dipole states would be shown 
in the forthcoming paper. 

For 10 Li, the present model successfully produces 
— 17.4 fm for the scattering length of the 2^ state 
in the 9 Li+n system as a signature of a virtual s- 
stateQ! HJ EH. The 1" state gives -5.6 fm, not so 
large negative value. Above the 9 Li+n threshold energy, 
twop-state resonances are obtained at 0.22 MeV and 0.64 
MeV for 1 + and 2 + states with the decay widths of 0.09 
MeV and 0.45 MeV, respectively. From these results, it 
is found that the Pauli-blocking naturally describes the 
inversion phenomenon of s- and p-states in 10 Li, in addi- 
tion to the 11 Li properties. 



IV. SUMMARY 

In summary, we have considered newly the tensor cor- 
relation in 11 Li based on the extended three-body model. 
We have found that the tensor and pairing correlations 
play important roles in 9 Li with different spatial char- 
acteristics, where the tensor correlation prefers a shrunk 
spatial extension. The tensor and pairing correlations in 
9 Li inside n Li are then Pauli-blocked by additional two 
neutrons, which makes the (Is) 2 and (Op) 2 configurations 
close to each other and hence activates the pairing inter- 
action to mix about equal amount of two configurations. 
As a result we naturally explain the breaking of magicity 
and the halo formation for 11 Li. We also reproduce the 
recent results of the Coulomb breakup strength and the 
charge radius of n Li. For 10 Li, the inversion phenomenon 
is explained from the Pauli-blocking effect. 

In this study, we focused on the tensor correlation, 
which is newly considered to figure out the s-p shell gap 
problem. However, the unified treatment of the effec- 
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tive interactions was not accomplished and is beyond 
the scope of this paper. This would require a consistent 
treatment of the short-range correlation in the realistic 
interaction while retaining the tensor force explicitly to 
describe the tensor correlation Ea. 
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